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A minor merger scenario for the ultraluminous X-ray 
source ESO 243-49 HLX-1 
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ABSTRACT 

The point-like X-ray source HLX-1 is the brightest known uhraluminous X-ray source 
and Hkcly the strongest intermediate-mass black hole candidate. HLX-1 is hosted 
in the SO galaxy ESO 243-49, but offset with respect to the nucleus, and its optical 
counterpart was identified with a massive star cluster. In this paper, we study, through 
A'^— body/smoothed particle hydrodynamics simulations, the scenario where ESO 243- 
49 is undergoing (or just underwent) a minor merger with a gas-rich low-mass late-type 
galaxy. The simulations suggest that the observed star formation rate (SFR) in ESO 
243-49 is a consequence of the interaction and that the companion galaxy already 
underwent the second pericentre passage. We propose that the counterpart of HLX-1 
coincides with the nucleus (and possibly with the nuclear star cluster) of the secondary 
galaxy. We estimate that, if the minor merger scenario is correct, the number density 
of X-ray sources similar to HLX-1 is w 10~^ Mpc^"^. 
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1 INTRODUCTION 

The point-like X-ray source 2XMM J011028. 1-460421 
(hereafter, HLX-1), with a maximum luminosity ~ 10*^ erg 
s"^ (Farrell et al. 2009, hereafter F09; Godet et al. 2009), 
is the brightest known ultraluminous X-ray source (ULX, 
see Feng & Soria 2011 for a review) and likely the strongest 
intermediate- mass black hole (IMBH) candidate (see van der 
Marel 2004 for a review). Estimates of the black hole (BH) 
mass based on the X-ray spectrum range from « 500 to 10* 
Mq (F09; Davis et al. 2011; Servillat et al. 2011). HLX-1 
is located in the outskirts of the SO galaxy ESO 243-49 (lu- 
minosity distance ~ 95 Mpc), ~ 0.8 kpc out of the plane 
and ~ 3.3 kpc away from the nucleus. HLX-1 has an opti- 
cal counterpart (Soria et al. 2010, 2012, hereafter SIO, S12, 
respectively; Farrell et al. 2012, hereafter F12), whose asso- 
ciation with ESO 243-49 is confirmed by the redshift of the 
observed Ha emission line (Wiersema et al. 2010). 

S«;i/f/Ultra-violet optical telescope (Swift/XJYOT) ob- 
servations of ESO 243-49 show asymmetric ultra-violet (UV) 
emission at ~ 2000 A (SIO, Webb et al. 2010). The UV 
emission is centred on the bulge of ESO 243-49, but has 
an asymmetric spatial distribution. A background galaxy at 
z ~ 0.03 can explain the asymmetry of the UV emission and 
contributes to a small fraction of the UV emission observed 
with Swift/UVOT (F12, S12). The UV emission centred on 



ESO 243-49 indicates ongoing star formation (SF) at a rate 
~ 0.03 M0 yr~^, consistent with a younger stellar popu- 
lation superimposed to the ~ 5-Gyr-old dominant stellar 
component of ESO 243-49 (SIO). This has been interpreted 
as the fingerprint of a recent minor merger between the SO 
galaxy and a gas-rich dwarf galaxy (SIO), whose nucleus may 
be identified with the counterpart of HLX-1. The presence 
of prominent dust lanes around the nucleus of ESO 243- 
49 (F12) is another hint of a recent gas-rich merger (e.g., 
Finkelman et al. 2010; Shabala et al. 2011). The possibility 
that minor mergers switch on hyperluminous X-ray sources 
(HLXs, i.e. ULXs with X-ray luminosity -assumed isotropic- 
Lx > 10*^ erg s~^) was initially proposed by King & Dehnen 
(2005, see also Bellovary et al. 2010). 

Furthermore, the optical counterpart of HLX-1 is asso- 
ciated with UV emission in both the near and the far UV, 
as shown by Hubble Space Telescope (HST) and Very Large 
Telescope (VLT) photometry (F12; S12). The UV emission 
from the counterpart of HLX-1 was also found to be variable 
(S12). 

The nature of the HLX-1 counterpart is still debated. 
Fits to the HST data (F12) indicate a total stellar mass of 
4 — 6 X 10^ M0, whereas both a very young (~ 10 Myr) and 
a very old (~ 13 Gyr) age are possible. An intermediate-age 
solution is not allowed. In the case of a very young stellar 
population minimal reprocessing from the disc is required, 
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whereas the old population explanation predicts that most of 
blue/UV emission comes from disc irradiation (F12). On the 
basis of VLT data, S12 show that the optical/UV component 
is variable and therefore dominated by disc irradiation, and 
exclude the scenario of a S> 10* M© young star cluster. The 
two remaining scenarios for the counterpart are a old ~ 10^ 
M© star cluster (likely a globular cluster) and a « 10* M0 
young star cluster. The latter explanation is complicated 
by the fact that such small star clusters unlikely host very 
massive (> 100 Mq) BHs (e.g., Portegies Zwart & McMillan 
2002). 

A further issue is represented by the observed X-ray 
variability of HLX-1, with a semi-regular period of ~ 380 
days (Servillat et al. 2011), and by the mechanism driving 
the accretion. Lasota et al. (2011, hereafter Lll) propose 
that the observed X-ray variability may reflect the orbital 
period of the companion star. In this case, the outbursts are 
triggered by periastron passages for a sufficiently eccentric 
orbit of the companion star (so that the tidal disruption 
radius is of the order of the periapsis distance). Considering 
the high mass transfer rate required to power HLX-1 (~ 
10~* Mq yr^^), Lll suggest that the companion star of the 
IMBH can be an asymptotic giant branch (AGB) star and 
that its mass cannot be too small (<< 3 M©). 

In this paper, we propose a new interpretation con- 
nected with the minor merger scenario (initially proposed by 
SIO). We investigate, through A'^— body /smoothed particle 
hydrodynamics (SPH) simulations, the possibility that the 
optical counterpart of HLX-1 is the nucleus of a gas-rich low- 
mass late-type galaxy, undergoing merger with ESO 243-49. 
In line with this interpretation, we suggest that HLX-1 is 
associated with the central BH of the disrupted galaxy and 
with its surrounding nuclear star cluster (NC). This entails 
that the counterpart of HLX-1 consists mainly of very old 
(> 10 Gyr) stars, but has also a younger population. In fact, 
NCs are known to have multiple stellar populations (Rossa 
et al. 2006; Walcher et al. 2006), generally dominated by 
a old component, but where a fraction of young stars is 
also present (e.g., the case of the Milky Way, see Pfuhl et 
al. 2011). Our simulations show that the overall scenario in 
which HLX-1 is the nucleus of a stripped, gas-rich low-mass 
galaxy is consistent with the available observations. 



1.1 The minor merger scenario 

A minor merger scenario between the SO and a gas-rich low- 
mass galaxy would naturally explain the UV emission and 
the dust lanes observed in the central region of ESO 243- 
49, as well as (at least) part of the UV emission in coin- 
cidence with the HLX-1 counterpart. Furthermore, a mi- 
nor merger scenario has other intriguing implications. For 
example, the counterpart of HLX-1 can be identified with 
the nucleus of the disrupting galaxy. This provides a pos- 
sible explanation for the mass of the BH in HLX-1, as the 
nuclei of dwarf galaxies have long been suspected to har- 
bour the low-mass tail of super-massive BHs (SMBHs, see, 
e.g., Barth et al. 2004; Peterson et al. 2005; van Wassenhove 
et al. 2010; Reines et al. 2011). Recent studies (Graham 
2012a, 2012b) indicate that, in low-mass galaxies, the mass 
of the central BH scales approximately with the square of 
the host spheroid mass (rather than linearly, as was usually 
thought, e.g. Marconi & Hunt 2003), implying that low-mass 



galaxies with a 10* — 10® M© bulge can host BHs with a 



10^ 



10' M0, consistent with the expected mass of 



the BH in HLX-1. Finally, nuclei of galaxies in this mass 
range frequently host NCs (e.g., Boker et al. 2002; Gra- 
ham & Guzman 2003; Cote et al. 2006; Graham & Spitler 
2009; see Boker 2010 for a recent review). The presence of 
a central BH and that of a NC do not seem to be mutu- 
ally exclusive: in addition to the Milky Way (e.g., Schodel, 
Merritt & Eckart 2009 and references there), Filippenko & 
Ho (2003) identified at least one galaxy (NGC 4395) host- 
ing both a SMBH and a NC, and Graham & Driver (2007) 
subsequently reported the existence of two additional such 
galaxies (NGC 3384 and NGC 7457). The sample of galaxies 
hosting both a SMBH and a NC was substantially increased 
by Gonzalez Delgado et al. (2008, 2009), Seth et al. (2008) 
and Graham & Spitler (2009). Unlike most of the Galactic 
globular clusters, NCs have often a complex SF history, with 
multiple episodes of SF (Rossa et al. 2006; Walcher et al. 
2006). This suggests the possibility that the colours of the 
HLX-1 counterpart are the result of multiple populations, 
dominated by old stars, but with a smaller contribution by 
very young stars (plus possibly emission from X-ray repro- 
cessing in the disc). 

We therefore simulate different scenarios for the merger 
between a SO galaxy and a gas-rich companion, to check 
whether the simulated kinematics and SF history match the 
observed properties of ESO 243-49 and of the counterpart 
of HLX-1. 



2 METHOD: A BODY SIMULATIONS 

The initial conditions for both the primary galaxy (i.e. the 
model of ESO 243-49) and the secondary galaxy are gener- 
ated by using an upgraded version of the code described in 
Widrow, Pym & Dubinski (2008; see also Kuijken & Dubin- 
ski 1995 and Widrow & Dubinski 2005). The code generates 
self-consistent disc-bulge-halo galaxy models, derived from 
explicit distribution functions for each component, that are 
very close to equilibrium. In particular, the halo is modelled 
as a Navarro, Frenk & White (1996, NEW) profile. We use 
an exponential disc model (Hernquist 1993), while the bulge 
is spherical and comes from a generalization of the Sersic law 
(Prugniel & Simien 1997; Widrow et al. 2008). 

Both the primary and the secondary galaxy have a stel- 
lar bulge and a stellar disc. The giant SO galaxy has no 
gas, whereas the secondary galaxy has an initial gas mass of 
1.38 X 10* M0, distributed according to an exponential disc. 
Therefore, the initial configuration of the secondary galaxy 
is consistent with a low-mass gas-rich disc galaxy. The total 
mass of the secondary is ~ 1/20 of the mass of the primary, 
classifying the outcome of the interaction as a minor merger. 
The masses of the various components and the characteris- 
tics lengths of the simulated galaxies are listed in Table 1. 

Here, we report the results of two different runs, in 
which the masses and scale lengths of each galaxy are the 
same (as described in Table 1), but their orbital proper- 
ties (impact parameter, relative velocity, orientation angles 
and total energy, listed in Table 2) are different. The main 
common orbital feature between the two runs is that the 
centre of mass (CM) of the secondary galaxy is assumed to 
lie approximately on the same plane as the disc of the pri- 
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Table 1. Initial conditions: masses and scale lengths. 



Table 2. Initial conditions: orbital parameters. 



Model properties 


Primary 


Secondary 


Orbital parameters 


Run A 


Run B 


DM Mass [IQii M©] 


7.0 


0.3 


Impact Parameter h [kpc] 


10.0 


10.2 


M*'' [IQio Mq] 


7.0 


0.2 


Relative velocity v^^i [km s^'^] 


200 


100 


/b/d 


0.25 


0.25 


e,^,i) [raAY 


7r/2, n, 


7r/2, 0, 2.94 


Gas Mass'' [10« Mq] 





1.38 


Initial distance D [kpc] 


200 


150 


Halo scale length"^ [kpc] 


6.0 


3.0 


Es [10^ km2 s-2] b 


0.38 


-1.65 


Disc scale length [kpc] 


3.7 


3.0 


Ls [lO'^ km s-l kpc] <= 


2.0 


1.0 


Disc scale height [kpc] 


0.37 


0.30 


ed 


1.003 


0.997 


Bulge scale length [kpc] 


0.6 


0.6 


Orbit « 


prograde 


retrograde 



^ M-t is the total stellar mass of the galaxy (including both 

bulge and disc), /b/d is the bulge-to-disc mass ratio. 

The primary has no gas, while the gas of the secondary is 

distributed according to an exponential disc, with the same 

parameters (scale length and height) as the stellar disc. '^ We 

name halo scale length the NFW scale radius Rs = i?20o/c, 

where -R200 is the virial radius of the halo (NFW 1996) and c 

the concentration (here we assume c = 12 for both galaxies). 



mary galaxy (in run A they initially lie exactly on the same 
plane, whereas in run B there is an initial shift of 2 kpc). 
This was required as the offset between the HLX-1 counter- 
part and the disc of the SO galaxy is relatively small (~ 0.8 
kpc) . We initially set the CM of the secondary galaxy at a 
distance of D = 200 kpc and 150 kpc from the CM of the 
primary in runs A and B, respectively (Table 2). These dis- 
tances are larger than two virial radii of the primary. The 
adopted orbits are nearly parabolic (run A is slightly hy- 
perbolic and run B is slightly bound) and with very high 
eccentricity (Table 2), in agreement with predictions from 
cosmological simulations (Khochfar & Burkert 2006). 

The particle mass in the primary galaxy is 2.5 x lO'' Mq 
and 5x lO'* Mq for dark matter (DM) and stars, respectively. 
The particle mass in the secondary galaxy is 2.5 x IC* Mq for 
DM and 5 x 10"' Mq for both stars and gas'^. The softening 
length is 0.1 kpc. We integrate the systems for 4 Gyr after 
the first pericentre passage. We simulate the evolution of the 
models with the A?^— body/SPH tree code gasoline (Wadsley, 
Quinn & Stadel 2004). Radiative cooling, SF and supernova 
(SN) blastwave feedback are enabled, as described in Stin- 
son et al. (2006, 2009, see also Katz 1992). The adopted 
parameters for SF and feedback are the same as used in 
recent cosmological simulations capable of forming realis- 
tic galaxies in a wide range of masses (e.g., Governato et 
al. 2010; Guedes et al. 2011), and in recent simulations of 
galaxy-galaxy collisions (Mapelli & Mayer 2012). 



3 RESULTS AND DISCUSSION 

We followed the interaction between the primary (SO) and 
secondary (gas-rich) galaxy for 4 Gyr after the first peri- 
centre passage (corresponding to ~ 5 Gyr starting from the 



^ We checked by running simulations with different resolution 
that our choice of particle masses does not affect significantly the 
dynamics of the system (e.g., by inducing appreciable spurious 
dynamical friction or ejections). We do not find significant differ- 
ences, apart from a change in the resolution of the SF. 



" For the definition oi 9, 4>, ip, see figure 1 of Hut & Bahcall 

(1983). In particular, 9 is the angle between the relative velocity 

vector Vjci and the symmetry axis of the primary disc, </> 

describes the orientation of Vj-ei projected in the plane of the 

primary disc and ^p describes the orientation of the initial 

distance vector D (between the CMs of the two galaxies) in the 

plane perpendicular to v^-^i. 
^ Es is the specific orbital energy , i.e. the total energy divided 
by the reduced mass fj, = mi m2/(mi -|- 012) (where mi and m2 

are the mass of the primary and of the secondary galaxy, 

respectively). Eg = — G M/D + v^^^/2, where M = mi + 1712 is 

the total mass of the two galaxies, G is the gravitational 

constant and D the initial distance between the CMs. 

'^ Lb is the modulus of the specific orbital angular momentum, 

i.e. the angular momentum divided by the reduced mass. 

•^ e is the eccentricity (e = [1 + 2EsL^/{G M)^]'^^^). 

'^ A orbit is classified as prograde/retrograde depending on the 

alignment/counter-alignment of the orbital angular momentum 

of the secondary galaxy with respect to the spin of the primary 

galaxy 



aforementioned initial conditions). The sinmlations indicate 
that the merger phase can be very long-lived (> 3 Gyr), if 
the orbital angular momentum (orbital velocity) of the sec- 
ondary is low (high), as the secondary spends most of its 
time out of the disc of the primary. We note that a system 
like ESO 243-49 plus the HLX-1 optical counterpart is not 
necessarily the result of the last stages of the merger, but 
can form almost anytime after the first pericentre passage'^, 
because of projection effects. The only limitation from the 
available spectroscopic measurements is that the relative ve- 
locity between the CMs of the two galaxies is not too large, 
as the measurement of the Ha line shows that the HLX- 
1 counterpart may be 'red-shifted' by « 170 km s"'^ with 
respect to the centre of ESO 243-49 (however, this measure- 
ment is quite uncertain, see Wiersema et al. 2010). 

This fact is evident from Figs. 1 and 2, that show the 
projected mass density of the overall stellar component of 
the two galaxies in run A (similar considerations can be done 
for run B). In Fig. 1, the density of stars has been projected 
in the plane that best matches the observed position of the 
HLX-1 counterpart with respect to ESO 243-49: the centre 
of the secondary is ~ 0.8 kpc out of the plane defined by 
the disc of the SO and ~ 3.3 kpc far from the nucleus of 
the primary. If we measure the relative line-of-sight velocity 
(wlos) between the CM of the secondary and of the primary 



^ The request that the two galaxies already had a first pericentre 
passage is necessary to justify the estimated SFR from the centre 



of ESO 243-49 (see Table 3). 
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Table 3. SFR and orbital properties in the simulations. 



Figure 1. Projected mass density of stars in run A at t = 2.6 
Gyr after the first pericentre passage. The two galaxies have been 
projected to the plane that best matches the observed projected 
position of the HLX-1 counterpart with respect to ESO 243-49. 
The scale is logarithmic, ranging from 2.23 Mq pc~^ to 2.23 X 10* 
Mq pc~^. The line-of-sight velocity of the secondary galaxy is 
'red-shifted' by ~ 230 km s~^ with respect to the CM of the 
primary. 




Figure 2. Projected mass density of stars in run A at i = 2.6 
Gyr after the first pericentre passage. The two galaxies have been 
projected so that the primary is seen face-on. The scale is loga- 
rithmic, ranging from 2.23 X 10~^ Mq pc'^ to 2.23 X 10* M 
pc ^. 



© 



galaxy according to this projection, we find i)los = 230 
km s~^ with respect to the CM of the primary, where the 
plus sign indicates that the secondary is receding from the 
observer. 

Fig. 2 shows that the apparent coincidence of the sec- 
ondary with the position of the primary when projected as 
in Fig. 1 is partially an effect of projection, as the centres of 
the two galaxies are ~ 60 kpc far from each other. Fig. 2 also 
shows the large bar (~ 13 kpc) in the primary, whose forma- 
tion was induced by the interaction and which produces the 
boxy shape of the bulge^. Finally, the tidal tails surrounding 





t 


SFR 


SFRscc 


d'i 


V-LOS 




[Gyr] 


[10-3MQyr-l] 


[lO-^Moyr-i] 


[kpc] 


[kms-l] 


Run A 


0.3 


1.0 


1.0 


89 


184 




0.5 


2.2 


2.2 


118 


128 




1.0 


1.9 


0.5 


154 


43 




1.5 


4.5 


0.4 


152 


-39 




2.0 


20.0 


0.1 


110 


-128 




2.5 


48.5 


0.2 


29 


306 




3.0 


22.1 


<0.1 


112 


73 




3.5 


12.0 


<0.1 


120 


-38 




4.0 


10.4 


0.0 


62 


-213 


Run B 


0.3 


1.1 


1.1 


66 


71 




0.5 


0.9 


0.8 


65 


-53 




1.0 


1.3 


<0.1 


54 


140 




1.5 


4.7 


<0.1 


26 


-278 




2.0 


8.1 


0.0 


71 


14 




2.5 


15.6 


0.0 


34 


296 




3.0 


12.9 


0.0 


87 


-12 




3.5 


6.6 


0.0 


7 


-370 




4.0 


5.2 


0.0 


72 


-34 



t: elapsed time since the first pericentre passage; SFR: total 

SFR; SFRscc: SFR in the secondary; dz: three-dimensional 

distance between the CMs of the two galaxies; fLOS^ relative 

velocity between the CMs of the two galaxies along the line of 

sight. A different line of sight is chosen for each snapshot, to 

match the observed position of the HLX- 1 counterpart with 

respect to ESO 243-49. 



the secondary are also apparent from Fig. 2. Figs. 1 and 2 
represent the state of the system at t = 2.6 Gyr after the 
first pericentre passage (and the first pericentre passage oc- 
curred ~ 700 Myr after the beginning of the simulation). At 
this time, the secondary already passed twice through the 
inner disc of the primary and it is receding again towards 
the apocentre. 

The last two columns of Table 3 show the three- 
dimensional (ds) distance and the relative line-of-sight ve- 
locity (i)LOs) between the CMs of the two galaxies at dif- 
ferent times i, assuming (for each snapshot) the projection 
that best matches the observed location of the counterpart 
of HLX-1. As defined before, the plus (minus) sign indicates 
that the secondary is receding (approaching) with respect 
to the observer (assuming that the primary galaxy is at 
rest with respect to the observer), i.e. it is red-shifted (blue- 
shifted). The modulus of wlos is almost always consistent 
with the observations (which suggest a red-shifting by ~ 170 
km s~^ with a large uncertainty, Wiersema et al. 2010), as a 
consequence of the assumed initial relative velocity between 
the two galaxies (Table 2). Therefore, the simulated systems 
spend most of their time in a kinematic state consistent with 
the available observations of HLX-1: more accurate spectro- 
scopic measurements of the HLX-1 counterpart are required 
to put stronger constraints. 



•^ A number of (currently unavailable) kinematic data (e.g., a 
rotation curve) are necessary to constrain the existence of a bar 



in ESO 243-49 (e.g.. Bureau & Athanassoula 1999; Athanassoula 
& Bureau 1999; Bureau & Athanassoula 2005). 
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Figure 3. Run A, top panel: three dimensional distance between 
the CMs of the two galaxies (da) as a function of time; bottom 
panel: SFR as a function of time for the entire simulation (open 
squares, blue on the web) and for the secondary galaxy (crosses, 
red on the web), i = is the time of the first pericentre passage. 
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Figure 4. The same as Fig. 3, but for run B. 

3.1 The star formation rate 

Further hints about the formation and evolution of the 
ESO 243-49 system can be derived from the study of the 
simulated SF rate (SFR). We remind that the UV observa- 
tions indicate a SFR~ 0.03 M© yr~^ centred in the bulge of 
the SO, although with an asymmetric distribution pointing 
towards the location of the HLX-1 counterpart (SIO). UV 
emission is also associated with the counterpart of HLX-1 



(F12; S12), but it may be totally or partially due to the 
reprocessing of the HLX-1 accretion disc (S12). 

The SF history from our simulations is shown in Figs. 3 
and 4, as well as in Table 3. In particular, the bottom panels 
of Figs. 3 and 4 show the SFR in runs A and B, respectively. 
The top panels of Figs. 3 and 4 show da as a function of time 
in runs A and B, respectively. It is important to compare the 
distance between the two galaxies with the SFR, to have 
an idea of the influence of the pericentre passages on the 
SF. In both runs A and B, there is no SF before the first 
pericentre passage (that is the absence of points at f < in 
Figs. 3 and 4 is due to the fact that no gas particles were 
converted to stars at i < 0). This indicates the robustness 
of our numerical SF recipes (see, e.g., Stinson et al. 2006) 
and the stability of our initial conditions, as no spurious SF 
is induced in the first stages of the simulations. SF starts a 
few Myr after the first pericentre passage. Initially, most of 
the global SF is concentrated in the gas-rich nucleus of the 
secondary. 

After --^1 — 2 Gyr (but this value depends on the or- 
bital properties), most of gas is tidally stripped from the 
secondary galaxy and accreted by the primary (see Fig. 5). 
Given the low angular momentum of the selected orbits, 
the gas fiows almost radially to the centre of the primary, 
where it starts corotating with the bar, forming a ring at a 
radius ~ 5 — 8 kpc, and is partially funnelled towards the 
centre of the SO. Here, a burst of SF takes place, delayed 
by ~ 1 — 2 Gyr with respect to the burst in the nucleus of 
the secondary. The SF burst in the centre of the primary 
is particularly strong after the second (or more) pericentre 
passage of the secondary. 

The SFR in the secondary galaxy reaches a few ^ 10~^ 
Mq yr^^ at its maximum (slightly after the first pericentre 
passage), whereas the SFR in the bulge of the primary can 
be as high as ~ 5 x 10~^ Mq yr~^ (Table 3), consistent with 
the Snnft/UYOT observations of ESO 243-49. The main dif- 
ference between the two runs consists in the fact that the 
SFR in the centre of the secondary is quenched more gently 
with time in run A (in which the velocity between the two 
CMs is higher, and there are just 2 pericentre passages in 4 
Gyr) , whereas it drops quite abruptly at i ~ 1 Gyr after the 
first approach in run B (in which the velocity between the 
two CMs is lower, and there are 5 pericentre passages in 4 
Gyr). 

Thus, the SF history derived from our simulations sug- 
gests that the ESO 243-49 system is in the late stages of 
a merger (after the second pericentre passage of the sec- 
ondary), in which most of the SF takes place in the bulge of 
the SO. On the other hand, if part of the blue and UV emis- 
sion from the counterpart of HLX-1 is due to young stars 
(rather than to the reprocessing in the accretion disc), the 
SF associated with the secondary must not be completely 
switched off. Therefore, the ESO 243-49 system might be in 
the stage between the first and the second pericentre pas- 
sage, when the secondary still retains a sufficient SFR. 

A scenario where the nucleus of the secondary still hosts 
a (relatively) young stellar population is preferred to explain 
some properties of HLX-1. In fact, the X-ray variability of 
HLX-1 (with a possible periodicity of 380 days, Godet et al. 
2009, Servillat et al. 2011) may be connected with a mod- 
ulated mass transfer due to tidal stripping of a star in an 
eccentric orbit around the massive BH (Lll). For this sce- 
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Table 4. Evolution of the stellar and gas mass bound to the secondary galaxy. 



t 

[Gyr] 


[106 Mq] 


[106 Mq] 


[W Mq] 


M-1 

[10'^ Mq] 


/\/rscc 
[105 Mq] 


[105 Mq] 


Run A 0.3 


4 


13 


15 


35 


1.1 


1.1 


0.5 


4 


12 


14 


34 


3.7 


4.7 


1.0 


2 


4 


14 


33 


1.9 


1.9 


1.5 


2 


4 


14 


33 


0.6 


0.6 


2.0 


2 


2.5 


14 


33 


0.4 


0.4 


2.5 


0.8 


1.4 


9.1 


22 


0.4 


0.4 


3.0 


1.2 


1.4 


9.1 


20 


0.1 


0.1 


3.5 


0.8 


0.9 


8.8 


20 


0.1 


<0.1 


4.0 


0.8 


0.9 


9.0 


19 


<0.1 


<0.1 


Run B 0.3 


4 


12 


12 


25 


2.3 


2.5 


0.5 


3.4 


11 


12 


25 


1.9 


2.0 


1.0 


0.9 


2 


7.2 


13 


1.5 


2.8 


1.5 


1.4 


2 


7.1 


13 


<0.1 


<0.1 


2.0 


0.1 


0.1 


3.9 


6.5 


0.0 


0.0 


2.5 


0.0 


0.0 


1.4 


3.2 


0.0 


0.0 


3.0 


0.0 


0.0 


1.4 


2.3 


0.0 


0.0 


3.5 


0.0 


0.0 


1.8 


1.8 


0.0 


0.0 


4.0 


0.0 


0.0 


0.4 


0.6 


0.0 


0.0 



t: elapsed time since the first pericentre passage; M!'^q gg and M|°^: total mass of gas in the inner 0.35 kpc and 1 kpc of the secondary 



galaxy, respectively; M^ 
A^ro.35 and M-- 



total mass of stars in the in the inner 0.35 kpc and 1 kpc of the secondary, respectively; 
total mass of young stars (< 200 Myr) in the inner 0.35 kpc and 1 kpc of the secondary, respectively. 



nario to reproduce the observed properties of HLX-1, the 
mass of the companion star cannot be too small (~ 3 Mq, 
Lll), although we stress that this requirement for the donor 
mass is the result of many assumptions (about the orbit 
of the star, the mass of the IMBH and the mechanism of 
accretion) . 



3.2 Evolution of the stellar and gas mass bound 
to the secondary galaxy 

We now focus on the stellar and gas mass that remains 
bound to the secondary galaxy at a given time (Table 4). 
M»°o.35 and Ml''\ are the total mass of stars within 0.35 kpc 
and within 1 kpc from the centre of the secondary galaxy, 
respectively, at a given time (Table 4). The mass within 1 
kpc gives an estimate of the stars that remain bound to the 
secondary galaxy, as the tidal radius is of the order of 1 kpc. 
Initially, Ml°\ is of the order of a few x 10* Mq . This is quite 
high with respect to the mass of ~ 10^ — lO'^ Mq , inferred for 
the optical counterpart of HLX-1. However, we stress that 
the values listed in Table 4 refer to the total stellar mass 
still bound to the stripped galaxy rather than to the (much 
denser) NC that can be hosted inside the nucleus: a NC 
cannot be dynamically resolved in our simulations (which 
have a softening of 0.1 kpc, i.e. approximately 100 times the 
characteristic size of a NC). 

The stellar mass within 0.35 kpc is reported for compar- 
ison with the VLT observations by S12 (which were taken 
with a seeing ^ 0.7 — 0.8 arcsec, corresponding to ^ 0.35 
kpc for the distance of ESO 243-49). We take the obser- 
vations by S12 as a reference, because the counterpart of 
HLX-1 is ~ 1 magnitude fainter in S12 than in the HST 
data by F12, indicating that the contribution of disc re- 



processing is lower. We converted the simulated values of 
M|°o 35 into an estimate of the apparent R, V and B band 
magnitudes in the Johnson system, by assuming that the 
luminosity distance of ESO 243-49 is 95 Mpc, by adopting 
a Chabrier IMF (Chabrier 2001), and by using the Padova 
GALaxies AnD Single StellaR Population ModELs4 (GAL- 
ADRIEL; Girardi et al. 2000) for a nearly solar metallicity 
{Z = 0.019) stellar population. We stress that this procedure 
applied to N-Body/SPH simulations cannot be very precise, 
as it is limited by the intrinsic resolution of SFR and stellar 
ages in the simulations. We find that, at i = 4 Gyr after the 
first pericentre passage, the apparent R, V and B magni- 
tudes associated to M^^g.ss are i? = 20 ± 1, V = 21 ± 1 and 
B = 22 ± 1 for run A^ and i? = 24 ± 1, V = 25 ± 1 and 
B = 26 ± 1 for run B. The corresponding values reported 
by S12 for the counterpart of HLX-1 are R = 24.71 ± 0.40, 
V = 24.79 ±0.34 and B = 25. 19 ±0.30. Thus, run B at t > 4 
Gyr after the first pericentre passage is consistent with the 
VLT observations of the HLX-1 counterpart. The strongest 
constraint comes from the R band, as the population associ- 
ated with the simulated secondary galaxy is predominantly 
old. We expect that part of the light in the B and V bands is 
contributed by disc reprocessing. Instead, the nucleus of the 
secondary galaxy in run A is still too bright at i = 4 Gyr to 
match the observed photometry of the HLX-1 counterpart. 
We conclude that either the orbit simulated in run A is not 
consistent with the evolution of HLX-1, or the elapsed time 
since the first pericentre passage is longer than 4 Gyr. 

Therefore, the comparison between the luminosity of 
the simulated galaxy and the observed VLT photometry 
gives important constraints for the minor-merger scenario. 
In particular, the properties of the HLX-1 counterpart can 
be reproduced only if (i) a time i > 4 Gyr elapsed since the 
first pericentre passage, or (ii) the orbit of the secondary 
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Figure 5. Projected mass density of gas in run A at t = 0.2 
Gyr (top panel) and at t = 2.6 Gyr (bottom panel) after the first 
pericentre passage. Tlie two galaxies have been projected so that 
the primary is seen face-on. The CM of the primary coincides 
with the centre of the frames, and it is marked by a white circle 
in the top panel. Each frame is 160 kpc per edge. The scale is 
logarithmic, ranging from 2.23 X 10~^ M© pc~^ to 2.23 X 10^ 
M0 pc-2. 



galaxy is at least as bound as the orbit adopted in run B 
(because a more bound orbit implies a faster tidal stripping, 
as the secondary galaxy undergoes more pericentre passages 
in a shorter time) , or (iii) the mass ratio between the primary 
and the secondary galaxy is initially > 20 (but galaxies sig- 
nificantly smaller than our simulated secondary galaxy are 
unlikely to host IMBHs at their centre). 

Furthermore, our simulations suggest that the counter- 
part of HLX-1 is surrounded by a diffuse stellar halo, as 
M»°i > A/*°o.35 for most snapshots. Thus, it is very impor- 
tant to assess whether the counterpart of HLX-1 is a com- 
pletely naked star cluster or is surrounded by such diffuse 
stellar halo. The existence of such halo is not immediately 
evident from either HST or VLT images (F12;S12). This may 
be an issue for our model. On the other hand, from our sim- 
ulations we expect the surrounding halo to be two (or more) 



magnitudes arcsec"'^ fainter than the central NC, i.e. very 
difficult to disentangle from the background contribution of 
the SO galaxy. We must be very careful in interpreting what 
occurs at the (barely resolved) very centre of the simulated 
secondary galaxy, as spurious numerical effects can affect 
density profiles and tidal stripping. In a forthcoming paper, 
we will investigate (through a wider grid of simulations) for 
which orbital parameters the merger scenario can reproduce 
the observed photometry of the HLX-1 counterpart. 

The total mass of gas in the inner kpc of the secondary 
galaxy after the first pericentre passage is ~ 1 — 2 x 10^ 
M0 (whereas most of the gas is already stripped in the tidal 
streams) and decreases by a factor of ^ 10 in the next 4 
Gyr (Table 4). Interestingly, the total mass of young stars 
bound to the secondary galaxy (defined as stellar particles 
younger than 200 Myr) is of the order of « 10^ M© in the 
first Gyr after the first pericentre passage and drops to « 10* 
M0 in the following Gyrs. We stress that most of the young 
stars bound to the secondary galaxy are located in the inner 
0.35 kpc, being much more concentrated than the old stars 
(Table 4). A young stellar component of ^ 10'* Mq is in 
agreement with the recent observations published by S12. 

3.3 Comparison with GalMer simulations 

The two high-resolution simulations described in the previ- 
ous sections provide important hints for ESO 249-43, but 
do not allow any statistically significant considerations. On 
the other hand, it is noticeable that two simulations with 
completely different orbital properties (Table 2) give similar 
results, indicating that the evolution we describe is quite 
common for interactions between gas-rich dwarfs and SO 
galaxies. 

Interestingly, the main features of our runs A and B (i.e. 
timescales, SF evolution and gas dynamics) are in agreement 
with the relatively low-resolution simulations of the GalMer 
database (Di Matteo et al. 2007; Chilingarian et al. 2010, 
and references therein). In particular, if we select out of 
the GalMer database the 36 runs describing an interaction 
between a giant SO galaxy and a gas-rich disc dwarf galaxy, 
we can make the following remarks, (i) More than half of 
the considered GalMer simulations (especially those with 
a retrograde orbit and the highest specific energy) develop 
the formation of a dense gas ring surrounding the bulge of 
the SO galaxy, very similar to the one showed in Fig. 5. (ii) 
The interaction triggers the SF, initially at the centre of the 
secondary galaxy and, at later epochs, even in the nucleus 
of the primary: the SFR reaches (on average) a value of 
0.08 — 0.1 Mq yr~^, higher than in our simulation. However, 
this discrepancy can be explained by the different resolution 
(mass resolution in our runs is a factor of ^ 10 higher), by 
the different SF recipes and especially by the fact that the 
total gas mass is a factor of ~ 4 higher in the GalMer sample 
than in our simulations, (iii) The average time spent by the 
simulated systems in a configuration similar to ESO 243-49 
(i.e. in a configuration where the nucleus of the secondary 
is at ^ 10 kpc projected distance from the centre of the 
SO galaxy and where SF occurs in the nucleus of the SO) is 
~ 700 Myr, with respect to an average total duration of the 
merger of ~ 3 Gyr (that is the average duration of mergers in 
the considered GalMer simulations). Thus, we can estimate 
that a SO galaxy interacting with a secondary galaxy is in 
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a 'ESO 243-49-like configuration' for a fraction /dose ~ 0.20 
of the duration of the merger. 

There are substantial differences between our simula- 
tions and the GalMer sample: the mass-ratio between the 
secondary galaxy and the SO galaxy is a factor of 2 larger 
in the GalMer database with respect to our simulations, the 
baryon-to-DM ratio is up to a factor of ~ 10 higher than 
in our simulations, the orbits span a slightly different pa- 
rameter space with respect to ours (in GalMer Es > 0, and 
Ls is always slightly higher than our values), the recipes for 
SF are quite different, and the mass resolution is a factor 
of ^ 10 lower. Taking into account all these differences, the 
agreement between the main features of our simulations and 
those of the GalMer database is noticeable. 



3.4 A statistical estimate 

In this Section, we estimate the density of ULXs similar to 
HLX-1 (hereafter, we will call them HLX-like sources, for 
simplicity), under the hypothesis that they are associated 
with a minor merger, and that they have an outburst X-ray 
luminosity > 10*^ erg s~^ (i.e., the minimum X-ray lumi- 
nosity for a source to be defined a HLX, see e.g. King &: 
Dehnen 2005). In particular, we assume that the host of 
a HLX-like source is a gas-rich low-mass galaxy, undergoing 
merger with a giant galaxy. The requests that the secondary 
galaxy is gas rich and undergoing merger are motivated by 
the fact that a gas-rich merger is needed to trigger SF. In 
addition, we assume that the secondary galaxy has initially 
a stellar bulge in the 10* — 10® M0 range, i.e. the mass range 
of bulges that might host a central BH with mass 10^ — 10^ 
Mq and a NC (Graham 2012a, 2012b). 

We define phost as the stellar mass density of bulges that 
host HLX-like sources. Phost can be estimated as 



Phost = 360 Mo Mpc 



-3 



/close \ //bH,NC\ // 



X 



0.71/ V 0.2 y V 1 y V 1 

Pb \ /'/morg^ f fMT\ ( /duty ^ 



/merg = 0.05, in agreement with available data and mod- 
els, although this number is uncertain (see, e.g., D'Onghia, 
Mapelli &: Moore 2008; Jogee et al. 2009, and references 
therein). 

/mt is the fraction of time that a BH similar to the one 
powering HLX-1 can spend in mass transfer with a com- 
panion star. There are no estimates of /mt for the case of 
HLX-1, and thus we adopt /mt = 0.02 from Blecha et al. 
(2006), which is the predicted time spent in mass transfer 
by a 100 — 500 M© BH hosted in a young star cluster. The 
actual value of /mt niay be very different from the assumed 
one, if the companion star is being tidally disrupted by the 
BH (Lll) and/or if the BH is in a different environment 
from a young star cluster, /duty is the duty cycle, i.e. the 
fraction of time that the BH spends in a HLX state (i.e., 
X-ray luminosity -assumed isotropic- higher than 10'*^ erg 
s~^). We assume /duty = 0.4 as, since its first detection, 
HLX-1 was observed above 10*^ erg s^^ for ~ 150 days dur- 
ing its ~ 380 day semi-regular period (Godet et al. 2009; 
Servillat et al. 2011). However, /duty is very uncertain, as 
we do not know the actual mechanism that powers HLX-1 
and the light curve of HLX-1 covers only two periods, which 
are not particularly regular. 

Therefore, the expected number density of HLX-like 



sources is 



riHLX ■^ 10"^ Mpc"''' 



0.03 y I 2.2 X 108 Mo Mpc- V V 0.05 J V 0.02 J I 0.4 J ' 



where h is the Hubble parameter (h — 0.71±0.025 according 
to the seven-year Wilkinson Microwave Anisotropy Probe 
data, Larson et al. 2011), /dose is the fraction of the merger 
timescale during which the nucleus of the secondary galaxy 
is sufficiently close to the disc of the primary (see previous 
section), pb is the mass density of bulges in the local Uni- 
verse (pb = 2.2 X 10* /iMq Mpc""', according to Driver et 
al. 2007). g is the mass fraction of bulges in the 10* — 10® 
Mo range (i.e., those that might host central BHs in the 
10^ - lO'^ M© range, Graham 2012a). We adopt g = 0.03, 
according to the Schechter formalism (Driver et al. 2007; 
Li & White 2009). /bh.nc is the fraction of bulges in the 
10* — W'' Mq range that host both the central massive BH 
and a NC. We adopt /bh.nc ~ 1, although this represents 
likely an upper limit (see, e.g., van Wassenhove et al. 2010, 
but Graham & Spitler 2009 indicate that /bh.nc may be 
of this order of magnitude). Finally, /gas is the fraction of 
galaxies with a bulge in the 10* — 10® M© range that host 
a significant amount of gas (~ 10* Mq). We put /gas = 1, 
as an upper limit, although this value is uncertain, /merg is 
the minor merger fraction in the local Universe. We assume 



h 

an 



phost 



36OM0Mpc-3 
3.7 X 10* Mq 



(Tlbulge) 



(2) 



where (mbuige) is the average mass of bulges in the 10* — 10® 
Mo range, according to the Schechter formalism (Li & White 
2009). This estimate suggests that there are « 4 HLX-like 
sources in a sphere with a 100-Mpc radius. Since we now 
observe one HLX-1, either we are missing a fraction of HLX- 
like systems, or (more likely) we are overestimating some 
of the parameters in equation (1). Likely, HLX-like sources 
have a duty cycle < 0.4 (e.g., because other sources are in 
different accretion regime with respect to HLX-1) and/or 
*MT < 0.02. Another possibility is that /bh.nc < 1 or that 
/gas < 1. Therefore, we can take our result as an upper limit, 
to be refined through new data. 



4 CONCLUSIONS 

HLX-1, hosted in the SO galaxy ESO 249-43, is the bright- 
est ULX known so far and the strongest IMBH candidate 
(F09). The optical counterpart of HLX-1 is a massive com- 
pact star cluster ( > lO** Mo, SIO; F12; S12), whose age is 
uncertain: /fST observations (F12) are consistent with both 
a very young (« 10 Myr) and a very old (> 10 Gyr) mas- 
sive star cluster. While the optical and UV variability of the 
counterpart seems to exclude a young star cluster with a 
mass higher than ~ 10'* Mq (S12), theoretical models sug- 
gest that the companion of the BH in HLX-1 cannot be too 
old, challenging even the old-cluster scenario (Lll). How- 
ever, we note that the available constraints on the mass of 
the donor strongly depend on a number of theoretical as- 
sumptions and that alternative accretion mechanisms still 
need to be investigated (e.g., instabilities in a radiation- 
pressure dominated disc, Lll). 
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In this paper, we studied, through A'^— body/SPH sim- 
ulations, the scenario where the SO galaxy ESO 243-49 is 
undergoing (or just underwent) a minor merger with a gas- 
rich low-mass disc galaxy. The simulations show that the 
observed UV emission and the corresponding SFR (~ 0.03 
M© yr~^) in the bulge of ESO 243-49 can be explained as 
a consequence of the interaction. From the comparison be- 
tween the observed SFR in the bulge of ESO 243-49 and the 
simulated SF history, we suggest that, if the UV emission 
from the optical counterpart of HLX-1 is mostly/entirely 
due to the reprocessing of the accretion disc, the ESO 243- 
49 system is currently at a stage after the second pericentre 
passage of the companion galaxy. In alternative, if part of 
the UV emission from the optical counterpart of HLX-1 is 
caused by a recent SF episode, the ESO 243-49 system might 
be in between the first and the second pericentre passage of 
the companion galaxy. 

We propose that the counterpart of HLX-1 coincides 
with the NC of the disrupting secondary galaxy. Recent 
studies (Graham & Spitler 2009; Graham 2012a) indicate 
that low-mass galaxies with a 10* — 10^ M© stellar bulge 
can host both a NC and a 10^ - 10^ Mq IMBH. The NC 
scenario explains many properties of HLX-1: NCs are often 
composed by multiple stellar populations, mostly old stars 
with a younger population superimposed. Under this hy- 
pothesis, the counterpart of HLX-1 consists mainly of very 
old (> 10 Gyr) stars, but the stellar companion of the IMBH 
can be a younger star. This satisfies the requirement that 
the mass of the companion star of the IMBH be not too 
small (<< 3 Mq) to reproduce the observed X-ray luminos- 
ity and variability (see Lll). To further check our scenario, 
we need new UV observations and a simultaneous monitor- 
ing of the X-ray and UV emission, to understand whether 
the UV emission from the HLX-1 counterpart is entirely due 
to disc reprocessing or is partially connected with a young 
stellar population. 

We compare the extrapolated luminosity of the simu- 
lated secondary galaxy with the observed magnitude of the 
HLX-1 counterpart (from the VLT observations reported by 
S12). We find that a very late merger stage ( ^ 4 Gyr from 
the first pericentre passage) , or a relatively bound orbit (ec- 
centricity ;$, 0.997), or a primary-to-secondary galaxy mass 
fraction ^ 20 are required to be consistent with the ob- 
served R magnitude of the HLX-1 counterpart. On the other 
hand, this constraint comes from many assumptions about 
the mass-to-light conversion in our simulations and may be 
affected by their spatial resolution. In a forthcoming paper, 
we will make a more accurate comparison with the observed 
photometry, by considering a wider grid of simulations (in- 
cluding a few higher resolution runs). 

If our scenario is correct, we expect the density of HLX- 
like sources to be ~ 10^'' sources Mpc^"^, corresponding to 
« 4 sources in a sphere of 100-Mpc radius. 

From our simulations, we find that the three- 
dimensional distance of HLX-1 from ESO 243-49 might be 
much larger (even by a factor of ~ 20) than the observed 
projected distance, if the merger is still ongoing. In fact, 
the existing spectroscopic measurements do not exclude a 
shift by ~ 170 km s~^, along the line of sight, between ESO 
243-49 and the HLX-1 counterpart. Therefore, new spectro- 
scopic measurements are required, to establish with better 
accuracy the relative velocity between ESO 243-49 and the 



HLX-1 counterpart, and to improve the constraints on the 
kinematics of the merger. In addition, there might be a num- 
ber of X-ray sources analogous to HLX-1 that have been 
missed or mis-classified because they are (in projection) far- 
ther out from the host galaxy. Therefore, it will be crucial 
to search for sources analogues to HLX-1, to test our and 
other possible scenarios. 
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